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ABSTRACT 

The  problem  considered  is  the  application  of  a  one -dimensional 
maximum  principle  to  second  order,  linear  differential  equations  of 
the  form 

u"  +  g(x)u'  +  h(x)u  =  f(x)  f or  a  <  x  <  b 
with  associated  general  boundary  conditions  to  obtain  functions  z  (x) 
and  z?(x)  such  that 

z2(x)  ^.  u(x)  ^  z^x) 
on  [a,b].   The  functions  f,g  and  h  are  assumed  to  be  bounded.   We  wish 
to  determine  the  behavior  of  the  solution  u(x)  on  [a,b]  and  also  to 
obtain  reliable  numerical  estimates  of  u. 

The  basic  concepts  in  the  theoretical  background  are  expanded 
versions  of  a  presentation  in  Protter  and  Weinberger  [Ref.  4]. 
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I.   ONE -DIMENSIONAL  MAXIMUM  PRINCIPLES 

If  u  is  a  real-valued  continuous  mapping  of  a  compact  set  X,  then 
u(X)  is  compact  and  by  the  Heine-Borel  Theorem,  u(X)  is  closed  and  bounded, 
Thus,  u  is  bounded.   Then  there  exist  points  p  and  q  in  X  such  that  u(q)  ^ 
u(x)  £  u(p)  for  all  x  in  X;  that  is,  u  attains  its  maximum  (at  p)  and 
its  minimum  (at  q) .   This  is  a  specific  example  of  a  general  maximum 
principle.   We  begin  our  discussion  with  a  study  of  some  specific  maxi- 
mum principles;  namely,  where  and  under  what  conditions  functions  of  a 
single  real  variable  which  satisfy  differential  inequalities  can  attain 
their  maximum  value.   The  approach  parallels  a  similar  investigation  in 

Protter  and  Weinberger  [Ref.  4]. 

2 
Definition:   For  u  =  (x)  in  C   [a,b],  define  the  linear  operator  L 

on  (a,b)  as  follows: 

LM  s  ^"  +  $(*)U' 

where  g(x)  is  any  bounded  function. 

2 
Lemma  1 :   If  u  ■  u(x)  in  C   [a,b]  satisfies  the  differential  in- 
equality L  [u]  )  0  on  (a,b),  then  u  cannot  attain  a  maximum  at  a  point 
interior  to  (a,b). 

Proof:  Assume,  to  the  contrary,  that  u  assumes  a  relative  maximum 
at  c,  a  <  c  <  b.   Then  u'(c)  -  0  and  u"(c)  ^  0>  hence 

LM    =  u"(c)+  4(c)v!(c)   =  H'(C)$0 
which  is  a  contradiction  of  L  [u]  }  0  on  (a,b).   Hence  u  can  attain  a 
maximum  only  at  the  end  points  a  or  b.  We  note  that  the  boundedness  of 

g  was  necessary  in  the  proof  of  Lemma  1  since,  if  lim  g(x)  ■  °°  ,  then 

ir*c 

L   [u]   would  be  undefined.   The  boundedness  condition  for  g  may  be 
JX-=C 

weakened.   The  fact  that  g(x)  is  bounded  on  every  closed  subinterval  of 


(a,b)  will  be  sufficient  for  our  purposes.   It  is  important  to  note  that 
g  bounded  on  every  closed  subinterval  of  (a,b)  does  not  imply  that  g  is 
bounded  on   [a,b].   Consider  g(x)  =  1/(1-1)  on  (-1,1),  and  g(x)  unbounded 

as  x  — *"  1. 

2 

Lemma  2:   If  v  ■  v(x)  in  C   [a,b]  satisfies  the  differential  in- 
equality L  [v]  K   0  on  (a,b),  then  v  cannot  attain  a  minimum  at  a  point 
interior  to  (a,b). 

Proof:   Note  that  -v  satisfies  the  hypotheses  of  Lemma  1  since 
differential  operators  are  linear;  that  is  L   [-v]  ■  -L   [v]  >    0. 
Hence  -v  cannot  attain  a  maximum  interior  to  (a,b)  which  implies  v  cannot 
have  a  minimum  interior  to  (a,b).   We  now  have  as  tools  two  simple  one- 
dimensional  minimum  and  maximum  principles.   For  convenience  we  will 
generally  discuss  explicitly  only  the  maximum  principles;  but  with 
reasoning  similar  to  the  proof  of  Lemma  2,  we  will  be  implicitly  dis- 
cussing the  related  minimum  principle. 

To  continue  our  development  we  relax  the  condition  of  strict  in- 
equality on  the  differential  inequality  and  now  consider  functions  u 
such  that  L   [u]    0.   The  solution  u  identically  a  constant  does  not 
satisfy  the  strict  inequality  in  the  hypotheses  of  Lemmas  1  and  2  since 
L   [u]  =  0.   We  now  allow  as  a  special  case  the  constant  function,  which 
attains  its  maximum  at  every  point  in   [a,b].   We  will  show  in  a  Corol- 
lary to  Theorem  1  that  the  functions  u  =  constant  are  the  only  func- 
tions which  satisfy  L   [u]^0  and  attain  their  maximum  interior  to  (a, 
b). 

In  investigating  the  behavior  of  a  function  u  which  attains  its 
maximum  at  x  ■  a,  we  might  imagine  that  u  could  have  a  graph  of  the  form 
shown  in  Fig.  la  (and  similarly  if  the  maximum  is  attained  at  x  =  b). 


That  is,  u  attains  its  maximum  at  x  -   a  with  u'(a)  =  0  (or  at  x  =  b 
with  u'(b)  =0). 


a. 


(a) 


b 


a. 


a 


(c) 


(b) 
Fig.  1 
We  show  now  that  this  behavior  is  impossible  if  u  satisfies  L   [u]  ^  0, 
If  M  is  the  maximum  of  u  on   [a,b],  then  u(a)  =  M  implies  u'(a)  <C  0 

(Fig.  lb):   and  u(b)  =  M  implies  u'(b)  >  0  (Fig.  lc) . 

2 
Theorem  1:   Let  u  ■  u(x)  in  C   [a,b]  be  a  nonconstant  function 

which  satisfies  the  differential  inequality 

V+  gOL)V  >,  0 

on  (a,b)  and  has  one-sided  derivatives  at  a  and  b,  and  suppose  further 
that  g  is  bounded  on  every  closed  subinterval  of  (a,b).   Then: 

(i)  the  maximum  of  u  occurring  at  x  =  a  and  g  bounded  below  at 
x  =  a  imply  u'(a)  K.     0; 

(ii)  the  maximum  of  u  occurring  at  x  =  b  and  g  bounded  above  at  x  =  b 
imply  u'(b)  >  0. 

Proof:   Suppose  that  u(a)  =  M  and  u(x)  ^  M  for  a  <  x  ^  b.   Since 
u  is  nonconstant  by  hypothesis,  there  exists  a  point  d  in  (a,b)  such 
that  u(d)  ^  M„   Define  the  function 

2(H)    ^    £        -   1 
with  A  a  positive  constant  to  be  determined.   We  note  that  z(x)  is 

positive  on  (a,b)  and  that 


LK1-  B*+$iti9'to  -  A  e         +  gooAe 


=  A  [  A  +  9^)J  e 


We  choose  A  large  enough  so  that  L  [z]  )  0  on  (a,d);  that  is,  we 
select  A  such  that  it  satisfies  the  inequality 

A  >  -  9(« 

on   [a,d]  which  is  possible  by  the  boundedness  assumptions  on  g(x). 
Now  form  the  function 

Wit)   ■  V-(TL)  +  Z  Ifti 
where  E  is  a  positive  constant  chosen  so  that  it  satisfies 

r    •        -ML-Wd) 

This  is  possible  since  the  numerator  is  positive  and  since  u(d)  \    M. 


Now 


LM  -  LM  +  ILM>o 


on  (a,d)  and  hence  the  maximum  of  w  on   [a,d]  must  occur  at  one  of  the 
end  points  by  Lemma  1.   We  have 

UJ(fi)   *  1A.OK)   4-  Z   Wa)     -  W. 

and 

upon  substituting  the  inequality  for  E;  hence  w(a)  X  w(d)  implies 
the  maximum  occurs  at  a.   Obviously,  the  one-sided  derivative  of  w  at 
x  ■  a  cannot  be  positive  for  then  w(a)  would  not  be  a  maximum.   Thus 


and  since  z ' (a)  =  A  y   0,  it  follows  that  u'(a)  \    0,  the  desired  re- 
sult in  (i). 

Suppose  that  u(b)  =  M  and  u(x)  ^  M  for  a  ^  x  ^  b.   Since  u  is 
nonconstant  by  hypothesis  there  exists  a  point  d  in  (a,b)  such  that 
u(d)  ^  M.   Define 

z(il)  =   e  -  1 

with  A  a  positive  constant  to  be  determined.   Then  z(x)  is  positive  on 
(a,b)  and 

LM  ■  2+  900  2  *  A  e         -  A  gwe 


-     .  ate-  y 

■  A  l  A  -  g«i    e 


We  choose  A  large  enough  so  that  L   [z]  )  0  on  (d,b);  that  is,  we 
select  A  such  that  it  satisfies  the  inequality 

A  >  $m 

on   [d,b]  which  is  possible  by  the  boundedness  assumptions  on  g(x). 
Now  we  proceed  as  before  and  form  the  function  w(x)  =  u(x)  +  E  z(x),and 
we  reach  the  conclusion  that  w   attains  a  maximum  at  x  «  b.   The 
one-sided  derivative  of  w  at  x  *  b  cannot  be  negative  for  then  w(b) 
would  not  be  a  maximum.   Then 

and  since  z'(b)  ■  -A  ^  0,  it  follows  that  u'(b)  >  0,  the  desired  re- 
sult in  (ii), 

Corollary  1A;   (One -Dimensional  Maximum  Principle) 

2 
Suppose  u  ■  u(x)  in  C   [asb]  satisfies  the  differential  inequality 
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u"  +  g(x)  u'  ^  0  on  (a,b),  with  g(x)  bounded  on  every  closed  subinterval 
of  (a,b).   If  u  attains  the  relative  maximum  M  at  an  interior  point  c 
of  (a,b),  then  u(x)  s  M. 

Proof:  Assume  u  has  a  maximum  at  c  in  (a,b),  then  u'(c)  ■  0.  Ap- 
plying Theorem  1  to  the  interval  (a,c)  and  (c,b),  we  conclude  that  u  is 
constant.   This  concludes  the  proof.   We  may  prove  Corollary  1A  in- 
dependently of  Theorem  1  as  follows: 

Assume  that  u(c)  ■  M  and  that  there  exists  a  point  d  in  (a,b)  such 
that  u(d)  <  M.   Take  d  >  c  for  convenience.   Then  define  the  function 

A        A 

where  A  is  a  positive  constant  to  be  determined.  We  note  that:  (i) 
z(x)  <  0  for  a  <  x  <  c,  (ii)  z(c)  -  0,  and  (iii)  z(x)  >  0  for  c  <  x 
<  b.   We  have 

A-Z  A- 1 


a[(a-i)  +  gwa-«)j  (t-*) 


We  choose  A  large  enough  so  that  L  [z]  s    0  for  a^x^  d;  i.e.,  we 
take  A  to  satisfy 

(A-i)  +  $M  (*-  a)  >  o 

or  equivalently 

A  >  -<$M(x--fi)  +  i 

which  is  possible  since  g  is  bounded.   Define  the  function 
where  E  is  a  constant  satisfying 
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°<E  <  -777T 


The  existence  of  E  is  guaranteed  by  the  hypothesis  on  the  point  d  and 
the  function  z.   Hence  on  (a,c) 

rw)  ~  nd)  +  JE  ?/jc)  <  A 

by    (i).      Also 

At    the   point  c 

Hence  w  attains  a  maximum  greater  than  or  equal  to  M  at  an  interior 
point  of  (a,d).   But 

LM  '  LM  t  £  Zf*3  >  o 

by  construction  of  the  function  z,  and  we  contradict  Lemma  1.   If 

d  <  c,  then  use  the  function 

A  A 

2(1)  «  (c-  a)     -   fltm) 

which  satisfies  (i)  z(x)  )  0  on  (a,c)  (ii)  z(c)  =  0,  and  (iii)  z(x) 

O  on  (c,b).   Note  that  the  inequalities  (i)  and  (iii)  are  reversed 
for  this  case.   This  concludes  the  proof.   It  would  also  be  possible  to 
use  the  function  M  .,,  -v 


as  the  auxiliary  function  in  the  proof  of  Theorem  1,  if  d  /  c,  and 
if  d  >  c. 


-AM-*)      A 


ii 


Corollary  IB;   A  nonconstant  u  satisfying  the  differential  in- 
equality L  [u]  ^  0  on  (a,b)  cannot  have  a  relative  maximum  at  an 
interior  point. 

Proof:   Assume,  to  the  contrary,  that  u  is  nonconstant  and  has  a 
relative  maximum  at  an  interior  point  c.   Then  apply  Corollary  1A  to  a 
subinterval  I  containing  c,  where  I  is  small  enough  that  u(c)  is  an 
absolute  maximum  on  I.   This  yields  the  immediate  contradiction  that 
u  S  M. 

Corollary  1C :  A  nonconstant  function  u  satisfying  the  differential 
inequality  L  [u]  ^  0  on  (a,b)  can  have  at  most  one  relative  minimum  in 
the  open  interval  (a,b). 

Proof:  We  note  that  between  any  minima  there  must  be  a  relative 
maximum  which  contradicts  Corollary  IB. 

The  boundedness  properties  of  g  are  required  for  Theorem  1  and 
Corollary  1A  as  shown  in  the  following  example.   Consider  the  functions 
u(x)  =  cos  x  and  g(x)  =  -cot  x,  which  satisfy  the  differential  equation 

and  here  u"  +  g(x)  u1  ^    0.   Conclusion  (i)  of  Theorem  1  does  not  hold  on 
[0,  IV  /2]  since  u'(0)  =  -sin(0)  ■  0  and  u  is  a  nonconstant  function. 
On  [-fK/2,  f^/2]  the  conclusion  of  Corollary  1A  fails  to  hold  since 
u(0)  =  cos(0)  =  1  is  a  maximum  at  an  interior  point  and  u  is  a  noncon- 
stant function.   We  note  that  g(x)  is  unbounded  in  the  neighborhood  of 
the  origin. 

Definition:   A  function  u(x)  has  a  horizontal  point  of  inflection 
at  x  =  c  if  u'(c)  =  0  while  u  is  strictly  increasing  or  strictly  de- 
creasing in  some  interval  containing  c. 
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Corollary  ID:  A  function  u  satisfying  the  differential  inequality 
L  [u]  ^  0  on  (a,b)  can  have  no  horizontal  point  on  inflection  in  (a,b). 

Proof:  Assume,  to  the  contrary,  that  at  the  point  c  in  (a,b),  u'(c) 
0.   Then  on  [a,c]  or  [c,b],  u  must  attain  its  maximum  on  the  interval  at 

x  =  c,  a  contradiction  to  Theorem  1. 

2 
Definition:   For  u  =  u(x)  in  C   [a,b]  and  h(x)  and  bounded  function, 

define  the  linear  operator  (L+h)  on  (a,b)  as  follows 

where  g  and  h  are  any  bounded  functions. 

■ 

In  this  more  general  setting,  it  is  necessary  to  modify  some  of  the 
preceeding  discussion  as  shown  by  the  following  examples.   Consider  the 
differential  equation 

%'V  -pt   s  0       on     [fijl^] 

which  has  the  solution  u(x)  =  sin  x.   Now  u  attains  a  maximum  of  1  at 
x  =  1^/2,  an  interior  point  of  [0,1^]. 

We  also  see  that  the  condition  h(x)  ^  0  is  not  sufficient  to  yield 
an  unrestricted  maximum  principle.   Consider 

#"-  n   *  0     cn    /-/>  (J 

X    ™x 

which  has  the  solution  u(x)  =  -e   -e   ,  and  u  attains  a  negative  maximum 

of  -2  at  x  =  0.   We  will  prove  that  a  nonconstant  function  u  satisfying 
(L+h)  [u]  ^  0  with  h(x)  ^  0  on  (a,b)  cannot  attain  a  non-negative  maxi- 
mum at  an  interior  point. 

Lemma  3 :   If  u  satisfies  (L+h)  [u]  £  0,  with  h(x)  ^  0  on  (a,b), 
then  a  nonconstant  function  u  cannot  attain  a  non-negative  maximum  at 
an  interior  point. 

Proof:   Assume,  to  the  contrary,  that  u(c)  =  M  for  a  ^    c  <  b  and 
u(x)  ^  M  on  (a,b).  Then  h(c)u(c)  ^=.0,  u'(c)  =  0,  u"(c)  ^  0;  hence 
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=  7t"fc)  +  hfc)  vte)  ±    o 

a  contradiction  of  (L+h)   [u]  >  0  on  (a,b). 

2 
Theorem  2 ;   If  u  =  u(x)  in  C   [a,b]  satisfies  the  differential 

inequality 

(Z+h)[*Q  *  % V gw n'+  /irxtit  ^O 

on  an  interval  (a,b)  with  h(x)  ^  0,  if  g  and  h  are  bounded  on  every 
closed  subinterval,  and  if  u  attains  a  non-negative  maximum  value  M  at 
an  interior  point  c,  then  u(x)  2  M. 

Proof:   Suppose  that  u(c)  =  M  and  u(x)  ^  M  for  a  <  x  <  b  and  M  ••*  0, 
Then  assume  that  there  exists  a  point  d  in  (a,b)  such  that  u(d)  ^  M  and 
for  convenience  take  d  ^  c.   Define  the  function 

where  A  is  a  positive  constant  to  be  determined.  We  note  that:  (i) 
z(x)  ^  0  on  (a,c):  (ii)  z(c)  ■  0;  and  (iii)  z(x)  )  0  on  (c,b).   Now 


(L+h)L*l  *    2*+  JM2  V  htt)  2 


M*-o)  At**)    .      f  Atx-c)     .  7 

Ae        +  Ajwe       +htt)^e       -  ij 


Ad-*) 


[a2+  Ayfr)+hm1  -  h(t) 
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We  choose  A  large  enough  so  that  (L+h)   [z]  )  0  on  (a,d),  that  is, 
so  that  A  satisfies  the  inequality 

r  "J       -Aft-*) 

Since  p      >0  and  h  ==  0  on  (a,b),  we  need  only  to  choose  A  such  that 

which  is  possible  since  g  and  h  are  bounded.  Now  form  the  function 

where  E  is  a  positive  constant  so  that 

Hence   on    (a,c) 

Tvtt)  *  Mxj  +  &  *wy   <  J/L 

since  z(x)  ^  0  on  (a,c).  And  by  the  definition  of  E, 

-?*////;  =  u(d)  +M-  KM)  ~  M 

Now 

iifid)  '  7f(C)  +  e*f<!)   m  M 

since  z(c)  =  0;  thus,  we  may  conclude  that  w  has  a  maximum  greater  than 
or  equal  to  M  which  is  attained  interior  to  (a,d).   But 

on  (a,d)  and  we  have  a  contradiction  to  Lemma  3  since  w(x)  is  noncon- 
stant. 

If  d  ^   c,  construct  the  function 

which  has  the  properties:   (i)  z(x)  ^  0  on  (a,c),  (ii)  z(c)  v   0,  (iii) 
z(x)  <^  0  on  (c,b).   Now 
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(L  +4)[b]  ±  *"  +  $ (i)  *'+  h/u  * 


Me  -  Ajfye         4  A/ty/e        -  lj 


«  e 


[A-toM+hm  -  ha) 


We  choose  A  large  enough  so  that  (L+h)  [z]  ^    0,  that  is  so  that  the 
following  inequality  is  satisfied: 

r  *  iii       Aa'c)  k  „ 

I  A-  Aim  +  h/x)  -  hm e  >  o 

since  £      y   0  on  (a,b).   And  since  h  ^  0,  choose  A  such  that 

A*-fyfc)+htt)  >o 

which  is  possible  since  g  and  h  are  bounded.   We  note  that  the  condition 
that  A  satisfies 

$K  A /$(*)/  -h  h/x)  >  6 

will  suffice  for  both  d  <*  c  and  d  ^>  c.   The  remainder  of  the  proof  is 
exactly  the  same  as  before  with  c  ^  d. 

Corollary  2A:   If  u  satisfies  (L+h)   [u]  £  0  on  (a,b)  with 
h  *£•  0  and  h(x)  3L   0,  and  u  assumes  a  non -negative  maximum  value  at  an 
interior  point,  then  u(x)  3E  0. 

Proof:   Let  M  ^  0  be  the  maximum  attained  by  u  on  [a,b].   By 
Theorem  2,  u(x)  £  M  and  hence 
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But  h(x)  £    0  and  h(x)  at  0  implies  M  m   0. 

2 
Theorem  3;   Suppose  that  u  =  u(x)  in  C   [a,b]  satisfies  the  dif- 
ferential inequality 

(L  +h)[\l  s  #"+  Saj  H'+h(*)K  *  0 

on  (a,b)  and  has  one-sided  derivatives  at  a  and  b,  that  h  ^  0  and  that 
g  and  h  are  bounded  on  every  closed  subinterval  of  (a,b).   Then: 

(i)   the  non-negative  maximum  of  u  occurring  at  x  ■  a  and  the  func- 
tion g(x)  4-  (x-a)h(x)  bounded  from  below  at  x  ■  a  imply  that  u'(a)  K.   0; 

(ii)   the  non-negative  maximum  of  u  occurring  at  x  =  b  and  the  func 
tion  g(x)  -  (b-x)h(x)  bounded  from  above  at  x  -  b  imply  that  u1 (b)  ^  0. 

Proof:  We  may  extend  the  proof  of  Theorem  1  by  noting  that  if  the 
maximum  of  uc  occurs  at  x  =  a,  then 

(l +h)[*      -l]  -  e      p  +  *$(*)  *  *>(*>  ft-  *      )] 


£  e        ltf+A$M+A  tie-*)  h/*)J 


Atx-«)r 

[' 

since  £  £  J'T  f°r  t  ^  0  and  h(x)  Z  0  by  assumption.   If  g(x)  + 
-t+l 

(x-a)h(x)  is  bounded  from  below  at  x  ■  a,  then  we  may  choose  A  such 
that 

If  the  maximum  of  u  occurs  at  x  ■  b,  then  as  above 


>. 


-  Att-k) 


[a  *--Ajm  +  a  fb  -x)  h(Jj 
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Hence  if  g(x)  -  (b-x)h(x)  is  bounded  above  at  x  "  b,  then  we  may  choose 
A  such  that 


A  >  d&J  -    (&  -  x)  hte) 


The  remainder  of  the  proof  is  analogous  to  Theorem  1. 

Corollary  3A:   If  u  j/  0  satisfies  (L4h)   [u]  ^  0  on  (a,b)  with 
h    0,  if  u  is  in  C2[a,b],  and  if  u(a)  £    0,  u(b)  ^  0,  then  u(x)  <  0 
on  (a,b). 

Proof:  Assume  that  there  exists  a  point  c  on  (a,b)  such  that  u(c) 
y     0.   Then  u(a)  £    0  and  u(b)  £*    0  imply  that  u  attains  a  non-nega- 
tive maximum  at  some  point  d  interior  to  (a,b),say  u(d)  =  M.  Then,  by 
Theorem  2,  u(x)  S    M  ^  0,  a  contradiction. 

A  useful  technique  in  determining  the  existence  of  a  maximum  (or 
minimum)  of  a  solution  to  a  differential  equation  at  an  interior  point, 
independent  of  the  preceeding  theory,  is  a  simple  direct  argument. 

Consider  the  problem 

which  does  not  have  the  format  of  the  problems  we  have  so  far  considered 

V. 

due  to  the  exponential  term  £  . 

We  note  that 


im*   k"+  e    *  -X  <  o 


on  (0,1).  Hence,  u  cannot  attain  a  minimum  in  (0,1);  otherwise,  if 
there  exists  a  point  c  in  (0,1)  such  that  u(c)  =  M  and  u(x)  ^  M  for 
x  ^r    c;  then  uM(c)  ^.  0  and  we  have 

on  (0,1),  a  contradiction. 

As  another  example  of  the  preceeding,  consider 

#" -2  6*2  (n')  -  i 
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whose  solution  u  cannot  attain  a  local  maximum  in  the  interior  of  any 
interval.   Let  I  be  an  interval  with  subinterval   [a,b].   If  there  exists 
c  in  (a,b)  such  that  u  attains  a  maximum  in   [a,b]  at  c,  then 

Hence 

0>  74n/d)  -  i+  2  aQ$[v'fc>"] }  ~  l+2&s/t>)  -  3; 

an  impossibility. 

We  may  even  apply  the  previous  results  to  boundary  value  problems. 
Consider 

for  Oxl,  with  u(0)  -  u(l)  ■  0.   Now 

U"+e\'<o 

on  (0,1);  hence  u  cannot  attain  a  minimum  interior  to  (0,1)  as  in  the 
first  example.   The  boundary  conditions  u(0)  ■  u(l)  =  0  imply  that  u'(0) 

.>  0  and  u'(l)  ^  0;  otherwise  u  would  attain  a  minimum  interior  to 
(0,1).   If  u'(0)  =  0  and  u(0)  ■  0,  then  from  the  differential  equation, 
u"(0)  ■  -1,   The  boundary  condition  u(l)  ■  0  forces  u  to  attain  a  mini- 
mum interior  to  (0,1),  a  contradiction.  We  obtain  a  similar  result  if 
we  assume  u(l)  =  u'(l)  ■  0,  since  u(0)  =  0.   Thus  u' (0)  >  0  and  u'(l) 

<  o. 
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II.   A  GENERALIZED  ONE -DIMENSIONAL  MAXIMUM  PRINCIPLE 

If  we  drop  the  restriction  that  the  function  h  must  be  non-positive 
we  may  then  obtain  a  maximum  principle  for  (L+h)  [u]  ^  0  by  introduc- 
ing a  function  w,  with  certain  desirable  properties,  such  that  the 

quotient  u/w  satisfies  the  previous  maximum  principles. 

2 
Theorem  4:   Suppose  u  =  u(x)  in  C   [a,b]  satisfies  the  differential 

inequality  (L+h)   [u]  ^  0  on  (a,b)  with  h(x)  bounded  and  g(x)  bounded 

below;  then  for  any  sufficiently  small  subinterval   [a',b']»  there  exists 

2 
a  function  w  in  C   [a,b]  such  that: 

(i)  w  >  0  on   [a',b'] 

(ii)   (L+h)   [w]  ^  0  on  (a',b'). 
Then  the  function  u/w  satisfies  the  maximum  principles  of  Theorems  2  and 
3  on  (a',b'). 

Proof:   Consider  the  function 

fr/x)  •  £  -  B(x-  «') 

for  a  ^  a'  (  x  ^  b'^  b,  where  B  is  a  positive  constant  to  be 
determined.  When  B  has  been  determined,  we  will  suppose  that  (b'-a1) 
is  such  that 

(i)  3(b'-*')3'<  L 

so  that  w    )0  on      [a',b'].      Now 

c  (.26)  +  JMh*6/*'*li]+  hit)  1-  B(x-a'f) 
-  -26 fd  /■  fr-a'JjM  +  {  (*-*')* hfr)l ]  +  kit) 
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Since  g  and  h  are  bounded  below,  choose  constants  G  and  H  such  that 
g  ^.  G  and  h  ^.  H.  We  now  make  a  restriction  on  (b',a')  so  that 

(2)  i  +fc~A?)Q  +  i  d-**f  H  >  o 

for  a'    ^    x    £    b*    and  choose  B   so   that    (L+h)    [w]    4s.    0  on   (a',b'), 
i.e.,  r  frlfij 

a  permissible  choice  since  h  is  bounded  above.  Hence  if  the  subinter- 
val  [a'jb'l  is  sufficiently  small,  the  function  w  with  properties  (i) 
and  (ii)  may  always  be  constructed. 

Define  the  function  v  -   u/w  on  [a',b'];  then 

=  (ri*"+  2Vu,>  twr'j-h^it)  /r'u>i  w'r)  +hMi~u> 


V"+[2  £  t  Wr'  +  i.  ti+^M  r  ±  0. 


Now  dividing  by  the  positive  function  w,  we  obtain 

Thus,  v  satisfies  the  differential  inequality 

on  (a',b')  where 

2 
By  the  construction  of  w,  w  in  C   [a',b'];  thence  g  and  h  are  bounded 

on  [a',b'].  Also  h  is  non-positive  by  conditions  (i)  and  (ii).   Thus 

v  satisfies  the  hypotheses  for  Theorems  2  and  3  on  subinterval  [a',b'], 

which  concludes  the  proof. 
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We  now  see  that  a  function  u  which  satisfies  (L+h)  [u]  ^  0  has 
certain  restrictions  on  its  zeros  on  the  interval  under  consideration. 

Corollary  4A;  If  the  function  u  satisfies  Theorem  4  on  (a',b') 
and  u(x)  re  0,  then  u  can  have  at  most  two  zeros  (between  which  u  is 
negative)  on  (a',b'). 

Proof:  We  may  take  the  interval   [a'^b*]  such  that  u(a')  =  u(b')  ■ 
0  and  u(x)  7P  0  for  a'  <  x  ^  b.  Assume  that  u(x)  >  0  for  a'  ^    x^  b. 
Then  v(x)  »  u(x)/w(x)  >  0  for  a'  <  x  ^  b '  and  v(a')  -  v(b')  m   0. 
Hence  v  has  a  non-negative  maximum  in  (a',b')  which  is  a  contradiction 
to  Theorem  4.  Thus,  u(x)  <^    0  on  (a',b').  Also,  there  could  not  be  more 
than  two  zeros;  otherwise,  if  there  exists  a  point  c  in  (a',b')  such  that 
u(c)  ■  0,  then  u(x)  <  0  on  (a',b')  forces  u'(c)  =  0.   But  then  u  must 
be  the  trivial  solution,  a  contradiction. 

Note  that  the  zeros  of  u  must  be  isolated;  i.e.,  if  u(a)  ■  0,  then 
u  cannot  vanish  on  some  finite  interval  to  the  right  of  a. 

Corollary  4B:   If  u  is  a  solution  to  the  differential  equation  (L+h) 
[u]  *  0,  then  u  can  have  at  most  one  zero  in  any  interval  (a',b')  where 
Theorem  4  holds. 

Proof:  Note  that  -u  also  satisfies  (L+h)  [-u]  ■  0.  By  Corollary 
4A,  u  and  -u  can  each  have  at  most  two  zeros  on  (a',b'),  between  which 
both  u  and  -u  must  both  be  negative,  an  impossibility.  Thus,  u  has  at 
most  one  zero  on  (a',b'). 

We  may  conclude  that  the  boundedness  hypothesis  in  Theorem  4  on  the 
function  g  is  essential  by  noting  the  following  example.   The  function 
u(x)  «  x  satisfies  u"  +g(x)u'  =  0  with  g(x)  ■  -2/x,  and  g  is  unbounded 
in  the  neighborhood  of  zero.   Consider  the  interval  of  the  form  (a,b) 
with  b  ■  0.   Then  u  attains  a  non-negative  maximum  of  zero  at  x  ■  b  and 
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u'(b)  =  0,  contradicting  conclusion  (ii)  of  Theorem  3,  hence  contra- 
dicting Theorem  4. 

Definition:   Let  r(x)  be  a  solution  of  the  differential  equation 

with  g  and  h  bounded  functions.   If  r  -jr  0  and  r(a)  =  0,  the  r  must  be 

T 

non-zero  on  some  finite  interval  to  the  right  of  a  by  Corollary  4B.   If 
r  has  another  zero  on  (a,  «^  ),  denote  the  first  such  zero  by  a*  and  call 
it  the  conjugate  point  of  a.   (If  r  JF  0  on  (a,  °°   ),  then  set  a*  =  o^>  ) . 

Thus  r  has  one  sign  on  (a, a*)  and  for  convenience  we  may  take  r(x) 
^   0  for  a  <^  x S  a*,  since  -r  satisfies  the  same  differential  equation. 
We  now  establish  a  necessary  and  sufficient  condition  for  the  existence 
of  the  function  w  of  Theorem  4  in  terms  of  a*. 

Corolllary  4C:  If  a*  is  the  conjugate  point  of  a,  there  exists  a 
function  w  ^  0  such  that  Theorem  4  holds  on  the  interval  [a,b]  if  and 
only  in  b  <^  a*. 

Proof:  Assume  that  b  ^    a*.   If  w  ^  0  on  [a, a*],  then  the  func- 
tion u/w  vanishes  at  a  and  a*  since  r(a)  -   r(a*)  =  0,  and  u/w  is  posi- 
tive on  (a, a*).   Hence  u/w  attains  a  non-negative  maximum  in  (a, a*)  and 
u/w  is  non-constant,  then  w  cannot  satisfy  (L+h)  [w]  ^  0  on  (a, a*)  with- 
out contradicting  Theorem  4.   Hence  by  a  contrapositive  argument,  we 
have  the  forward  implication  that  the  existence  of  an  acceptable  w(x) 
implies  b  <^  a*. 

To  get  the  reverse  implication,  let  b  in  (a, a*).   Note  that  r(x)  is 
bounded  from  below  by  a  positive  number  m  ■  m(c)  on  any  subinterval  [c,b] 
with  a  <  c.   Let  &(%—&) 

*m  *  2-  e 

where  A  is  a  positive  constant  to  be  determined.   Consider  the  function 
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ivii)*  W/J  +  B  *M 

Now  choose  E  small  enough  so  that  w  )  0  on  [a,b].   Now 

=  Ej(-Ae       )+y*)(-At      )  +  h  (*■)  (2-t       Jj 

We  want  to  choose  A  such  that  (L+h)  [w]  ^  0  on  [a,b]  which  can  be  done 

as  follows: 

At* -a)  AM- *) 

Now  -Ee      ^  0  since  E  and  e        are  positive.  We  make 
two  restrictions  on  A: 

-At*-*)  v  A    v     /n  Z 

(i)  1  -  2e        >  0,  i.e.,      ft     ' 


Since  h  is  bounded,  we  may  also  choose  A  large  enough  so  that 

Then  we  have  (L+h)  [w]  4s>  0,  and  w  is  a  function  for  which  Theorem  4  holds, 

If  h(x)  is  unbounded  or  g(x)  is  not  bounded  below  there  may  be  no 
interval  [a'jb1]  for  which  Theorem  4  holds.   Consider  the  function 

which  satisfies 

71  \  xT¥7t    =  C 

Now  as  x — ►  0+,  h(x)  =  x — *•  +  ****.   There  can  be  no  function  w  ^  0 
with  the  property  that  u/w  satisfies  the  maximum  principles  in  any  of 
the  intervals  [0,  l/n'fK  ]  (n  =  1,2,...)  since  u(l/n4T)  -  0  and  1/nlr-K) 
as  n— *•  c*-°    .   Thus,  for  all  £  >  0,  the  interval  (0,  £  )  contains  a 
zero  of  u  and  Corollary  4A  denies  the  existence  of  an  appropriate  w(x). 
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As  in  the  examples  in  the  previous  section,  we  use  a  direct  argument 
to  conclude  that  no  solution  of 

can  attain  a  minimum  value  greater  than  1  or  a  maximum  value  less  than  1. 
Assume  there  exists  a  point  c  such  that: 

(i)  u(c)  =  maximum;  hence  u"(c)  j£  0  and  we  have 

MM) 

e      —    e    ■  -  iiu(e)   ^  o 

which  implies  e   >e,  and  we  have  the  maximum  value  of  u  is  u(c)  ^  ln(e) 

-  i; 

(ii)  u(c)  =  minimum;  hence  u"(c)  £.   0  and  we  have 

e     -   e   ■  u  (c)  £=  0 

which  implies  that  e      «£•  e  and  we  have  the  minimum  value  of  u  is 
u(c)  £    ln(e)  -  1. 
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III.   INITIAL  VALUE  PROBLEM 
Consider  the  solution  to  the  differential  equation 

CD      #*+gflOJt/+/)/*in  --fte) 

which  satisfies  the  initial  conditions 

(2)        lite)*  %     )      W'M  =  % 

where  the  functions  f,  g  and  h  are  defined  on  the  interval  (a,b)  with 
g  and  h  bounded  and  "     and  0      given  constants. 

The  existence  of  a  solution  follows  from  the  classical  theory  of 
ordinary  differential  equations  [Ref.  3].  We  may  prove  the  uniqueness 
of  such  a  solution  independent  of  the  classical  theory  using  Theorem  4. 

Theorem  5:   Suppose  the  u  (x)  and  u  (x)  are  solutions  of  (1)  which 
satisfy  the  same  initial  conditions  (2)  in  an  interval  (a,b).   Then 
ul   S  u2  in  (a,b). 

Proof:  Define  u(x)  »  u. (x)  -u_(x).  Then  u  satisfies 

with  initial  conditions 

It  (A)*   U'fr)  ~0 
Now  assume  that  u  p  0  in  (a,b).   Then  by  Theorem  4  there  exists  a 
positive  function  w  on  an  interval  (a,a+£  )  such  that  the  maximum  of 
u/w  in  (a,a+£  )  must  occur  at  one  of  the  endpoints.   We  note  that 
is  a  positive  constant  depending  only  on  the  bounds  of  g  and  h.   I.e., 
(a+  £  )  -  £     *  £.  must  satisfy  conditions  (1)  and  (2)  in  the  proof  of 
Theorem  4.   But  -u  also  satisfies  the  same  equation  with  the  same  initial 
conditions,  hence,  also  by  Theorem  4,  -u/w  must  also  attain  its  maximum 
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at  one  of  the  endpoints,  a  or  a+.   Hence  u/w  must  attain  either  a  maxi- 
mum or  a  minimum  at  a.   But  at  x  ■  a 


/-u.\  I  -k'u/-    WW1 


since  u(a)  ■  u'(a)  =  0.   Now  the  function  u/w  satisfies  Theorem  3;  hence 
u/w  must  be  identically  a  constant,  and  from  the  initial  condition,  u=0 
on  [a,a+*2  ].   In  particular 

t4ia+€)  *  zt'fa+e)   -  £ 

Hence,  by  the  analogous  argument,  u  Er  0  on  [a+£   ,  a+2£.  ]  with  £   un- 
changed since  the  choice  of   £   is  independent  of  u.   Now  (b-a)  <  <**>   ; 
hence,  after  a  finite  number  of  steps,  we  may  conclude  that  u  35  0  on 
[a,b].   Thus  u.  &.    u~. 

Consider  the  problem 

?,'/_  J.     U'    =   £> 
with  the  initial  conditions 

2VtfJ  =    it '/»)  =  £> 

We  note  that  U[  =  0  satisfies  the  equations  trivially  and  for  ^   =  x2, 

<-{*'=    2-  i'(z%)    =0. 

But  u  ^t  u  on  (0,1).   Since  g(x) — >■    '=**       as  x — >■   0,  we  may  conclude 
that  the  boundedness  of  g  is  necessary  in  Theorem  5. 

We  can  easily  find  a  function  h(x)  such  that  the  equation 

it**  km  *  -  o 

has  two  solutions  satisfying 

(4/0)   =r  (I  '(0)   s?    O 

Obviously  u,  "S.   0  is  a  solution.   We  must  then  pick  h(x)  unbounded  in 

the  neighborhood  of  the  origin;  otherwise  Theorem  5  would  imply  that  the 

2 
trivial  solution  was  the  only  solution.  Choose  h(x)  =  -6/x  ,  then  u(x) 

3 
■  x  is  an  appropriate  solution  satisfying  the  initial  conditions. 
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Hence  we  see  that  the  boundedness  hypothesis  on  h  is  also  necessary  in 
Theorem  5. 

An  interesting  application  of  the  Mean  Value  Theorem  arises  in 
the  initial  value  problem 

(3)  w"+  e    -  1        for    x>o 

satisfying 

(4)  7//£j=  /  and    -y'/o)  =  °- 

Assume  that  u.  and  u.  are  solutions  of  the  differential  equation  (3) 
satisfying  the  same  initial  conditions  (4).  Then  the  function 

Ttu)  =  %  m  -  74  ft) 

satisfies  the  homogeneous  differential  equation 
By  use  of  the  Mean  Value  Theorem 

e  -  e     =  C%~uz)e  ;  o<e<\ 

and  the  differential  equation  becomes 
with 

(1(0)  *    7f7c)  ^O 

Since  h(x)  *  -e  fr  0  and  is  bounded  on  any  finite  sub- 

interval  of  (0,  o&  ),  we  may  apply  Theorem  5.  The  trivial  solution 
u  E  0  satisfies  the  above  initial  value  problem,  thus  u  is  unique 
and  u  £  u  .   Hence  we  have  the  solution  of  (3)  satisfying  (4). 
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IV.   BOUNDARY  VALUE  PROBLEM 

We  will  first  be  concerned  with  boundary  value  problems  of  the 
type 

(i)     rt'+gMu'+hM-K  =  fa) 

on  the  interval  (a,b)  subject  to  the  boundary  conditions 

(2)     *(*)-  %    ',  ka/-  yt  . 

We  may  derive  a  uniqueness  theorem  by  the  use  of  Theorem  4  for  solutions 
of  (1)  which  satisfy  (2).   The  sign  of  the  function  h(x)  has  great  signi' 
ficance  as  seen  in  the  simple  example: 

u"  +  n  =  o 

x-/hich  has  the  solutions  u  «*»  0  and  u  =  sin  x  on  the  interval  [0,^  ]. 
Both  u  and  u  satisfy  the  boundary  conditions  u(0)  =  u('IN')  =  0. 

Theorem  6:   Suppose  that  u  (x)  and  u  (x)  are  solutions  of  (1) 
which  satisfy  the  boundary  conditions  (2),   If  h(x)  ^  0  on  (a,b),  then 

Ul  ~  V 

Proof:      Set  u(x)    ■  u    (x)    -  u    (x) ;    then  u   satisfies 

■H"+4W  U'+  h/JL)  l{    —  O 

with  the  boundary  conditions 

K{<£)  -  ute)  =■  6 

By  Corollary  2A*  u(x)  -  0  on  (a,b).   By  analogous  reasoning  -u  satis- 
fies the  above  equations;  hence  -u  —  0  on  (a,b).   Thus  u  3t  0  on 
(a,b). 

In  order  to  incorporate  a  larger  class  of  boundary  value  problems 
we  consider  more  general  boundary  conditions  which  have  (2)  as  a  special 
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case.  We  consider  solutions  of  (1)  which  satisfy  the  boundary  condi- 
tions 

r—v'{*)00S6    +     l//*)^n&    -  % 

C3)  \  , 

I  n'fbJGost  +   Wh)  dm  4    ■  € 

where  £>-  ,  <^  ,  ^  .  and  '/*_  are  prescribed  constants  with  0  £&£3£    , 
0  £  4  ±   0-        .We  have  (2)  when&=  ^   ■'%  * 

To  derive  the  general  boundary  conditions  (3)  consider  the  boundary 
conditions 

/},  (-*'(*))    +   B,  y/aj  -  £j 

where  A  ,  B.,  C  ,  (i  ■  1,2)  are  given  constants.  We  may  normalize  the 
above,  i.e., 

y  am?  (  iWi  e?  Va*+*,2         i 

and  ^        ' 

Then  we  define  the  angles  &   and  <p   so  that 
and 

*****  Wv    '  'ty'  tffiv 

Then  0  5^^  and  0  £^i  TK   • 

Theorem  7 ;   Suppose  that  u  (x)  and  u_(x)  are  solutions  of  (1) 
which  satisfy  boundary  conditions  (3).   If  h(x)  —  0  on  (a,b),  then 
u.  £=  u-  unless  h  =r-  0,  6?"  =  ^  =  0;  in  which  case  u1  and  u_  may  differ 
by  a  constant. 

Proof:  As  before  set  u  ■  u  -  u_,  so  that  u  satisfies 
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and 

(4)    -K'(*)  &S6>    +   MM)  $vn  &  -    0 

(5>   7/W  out 4  f  n(k>i  $m4>  =  o 

Consider  u£  H,  a  non-zero  constant.   Then  (4)  and  (5)  imply 

Now  M  ■£  0;  hence  sin£^  =  sin^  ■  0  and  &  ■  <B  =  0.   By  Corollary  2A, 
M  ^  0  if  and  only  if  h  ss  0.   Thus,  if  u  and  u  differ  by  a  constant, 
then  h  S  0  and  Q  -  <£    -   0. 

Conversely,  assumed  =  <f>    ■  0  and  h  *£  0;  then  u  sE  M  satisfies  (4) 
and  (5)  trivially  and 

Hence,  u  and  u.  differ  by  a  constant. 

Next  assume  u  is  a  non-constant  function  which  is  positive  at  some 
point  in  [a,b].   By  Theorem  2,  u  attains  its  positive  maximum  at  a  or 
b.   Suppose  the  maximum  occurs  at  x  *  ft.   Then  by  Theorem  3,  u'(a)  ^  0. 
Now  u(a)  >  0and0^9^^',0^<^^J'   ;  hence 

-  //  '{*)  &jo  -h  -y/a)  £**,  e  >  o, 

since  cos  ^  and  sin  £  cannot  vanish  simultaneously  and  are  positive 
on  the  interval  under  consideration.   This  inequality  is  a  contra- 
diction of  (4).   If  the  maximum  occurs  at  x  =  b,  we  have  u'(b)  ^    0, 
u(b)  >  0,  and 

?■(  Yi)  tos  tf  -h  u/h )  £m  <fi  >  o 

a  contradiction  of  (5).   Hence  u  is  a  non-constant  solution  which  can 
never  be  positive.  Analogously  -u  is  also  a  solution  which  satisfies 
(4)  and  (5);  and  we  may  conclude  that  u  can  never  be  negative.  Thus  u 
2?  0  on  [ a , b  j . 
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Now  let  the  only  assumption  on  the  function  h(x)  be  that  it  is  bound- 
ed on  the  interval  under  consideration;  i.e.,  we  remove  the  restriction 
that  h(x)  be  non-positive.  We  will  now  prove  a  uniqueness  theorem,  but 
as  the  following  example  shows,  the  interval  considered  is  of  vital  im- 
portance.  The  simple  equation 

p*+  7i   «  O 

with  u(a)  ■  u(b)  =  0  has  only  the  trivial  solution  for  b-a  <  1"K  ;  and 
for  b-a  ■  IK  uniqueness  failed  as  was  shown  at  the  beginning  of  this 
section. 

Theorem  8 :   Suppose  that  u  (x)  and  u.(x)  are  solutions  of  (1) 
which  satisfy  the  same  boundary  conditions  (2)  with  h(x)  £=.    0  on  [a,bj. 
If  b  ^    a*,  where  a*  is  the  conjugate  point  of  a,  then  u.  S  u  . 

Proof:   Let  u  =  u-  -  u_;  then  u  satisfies  the  equation 

with 

Since  b  ^  a*  by  hypothesis,  Corollary  4C  implies  that  there  exists  a 
function  w  ^  0  on  [a,b].  Hence  by  Theorem  4  the  function  v  ■  u/w  satis- 
fies the  maximum  principles  of  Theorems  2  and  3  and  their  Corollaries, 
Now  v(a)  ■  v(b)  =  0  since  u(a)  ■  u(b)  =  0,  hence  v(x)  £  0  on  [a,b]; 
otherwise,  v  would  attain  a  positive  maximum  in  (a,b),  a  contradiction 
to  Theorem  2  and  Corollary  2A.  Since  w  is  a  positive  function,  u(x)  <£  0 
on  [a,b].  We  may  apply  the  same  argument  to  -u  since  it  satisfies  the 
same  differential  equation  and  boundary  conditions,  and  we  conclude  that 
-u  &    0.   Thus  u  S  0. 
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V.   APPROXIMATION  IN  BOUNDARY  VALUE  PROBLEMS 

We  are  concerned  with  the  same  boundary  value  problem  as  in  Section 
IV.   Thus,  we  seek  a  solution  of  the  differential  equation 

(i)  (M)fii]  2  2/V^//j^'+  hlKiu  *  /#; 

on  the  interval  (a,b)  which  satisfies  the  boundary  conditions 

(2)      7im  -  %    ;    W  (h)  •  "fj.  . 

Even  in  many  of  the  simplest  problems  an  explicit  expression  of 
the  solution  is  impossible  to  determine.  We  must  then  resort  to  series 
solution  techniques  or  other  essentially  numerical  methods.  Obviously 
an  approximation  solution  technique  which  yields  an  explicit  bound  for 
the  error  compared  to  the  exact  solution  is  desirable,  It  is  possible 
to  find  such  a  bound  by  finding  explicit  approximation  functions  which 
provide  an  upper  and  lower  bound  for  the  actual  values  of  the  solution. 

In  this  section  we  will  develop  necessary  conditions  for  the  exis- 
tence of  such  approximation  functions,  and  in  Section  VI  we  will  in- 
vestigate more  closely  the  approximation  functions  and  related  results. 

We  assume,  as  before,  that  the  functions  f,  g  and  h  are  bounded  on 
the  interval  under  consideration.   With  these  assumptions  it  is  possible 
to  determine  bounds  for  a  solution  u  of  (1)  and  (2)  using  the  maximum 
principle  in  Theorem  3  without  any  further  knowledge  about  u. 

Theorem  9;   Suppose  that  u  is  a  solution  of  (1)  satisfying  boundary 
conditions  if/.),  and  tnat  h(x)    0.   If  the  function  z  (x)  satisfies: 


(3)    (L+b)[iJ±4tl)      on     fab) 
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and  the  function  z-(x)  satisfies: 

(5)  (A+h)C*J  k  J-/*)      <W7   (tik) 

(6)  2t/«j*r,    wr/    V<y/& 

then 

for  a    fr    x  £    b. 

Proof:  Consider  the  function 


Then  v1  satisfies 


on   (a,b)  with 

V}fr)    &Q      tfyid     VJ(b)  ±  C. 

Apply  Corollary   3A;   v. (x)   £  0  on   [a,b].      Thus, 

V/t)  <r  *,  ft)  . 

Similarly  the  function 

-n  it)  ?  2$/&  -  *(t) 

satisfies 

on  (a,b)  with 

Thus  v-(x)  ^   0  on  [a,b]  and  we  have 

on  [a,b]. 

A  more  general  two-point  boundary  value  problem  is 

c?)    (JLtf)F*3-  W+ 3fcjit'+ ka)  u.   --ftic) 

on  (a,b)  where  the  solution  u  satisfies 

It '(b)  <?os0    +  u/t>)   St*\<fi  -   t-L 
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The  quantities  O    and  <p  are  pre-assigned  constants  wuch  that 

O^  &  ^  lr/2.      and      o  ±  <fi  ^  %   . 

Continuing  to  assume  h(x)  ^-  0,  we  may  extend  the  proof  of  Theorem  9 
to  get  the  following. 

Theorem  10;   Suppose  u  is  a  solution  of  (7)  satisfying  boundary  condi- 
tions (8).   Suppose  also  that  h(x)  ^  0  and  &■    ,  «p  are  in  [0,^/2] 
and  that  not  all  the  equalities  h  3  0,  ©  =  0  ,<p  =0  hold.   If  the  func- 
tion z  (x)  satisfies: 

(9)     (h+h)[$D  -  J't)       $h      f*jb) 
(10)        2;(bJ&scf>    +     */b)  **t4    ±    %L 

and  the  function  z_(x)  satisfies: 
(12) 

*±'/6)  (to/ 4  +  &W  ^4  *r?L 

then 

^m  ±  wit)  ±  tyx) 

Proof:   Set  the  function  v.  ■£  u  -  z  .   Then  v  satisfies 
and 

r/Wyutf  f  ft(b)  &o*4  -o  \ 

If  v.  is  always  positive,  then  Theorem  2  states  that  the  positive  maxi- 
mum of  v  occurs  at  either  a  or  b.   If  it  occurs  at  a,  then  v,(a)  ^   0 
and  v.  (a)  >  0.    But 
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can  occur  only  if  sin  Q    =  0,  i.e.,  ^  ■  0  and  v'(a)  ■  0.   Then  Theorem 
3  states  that  v..  (x)  is  identically  a  positive  constant.  Thus  h(x)  &  0 
by  Corollary  2A.  Similarly  v  cannot  attain  a  positive  maximum  at  b 
unless  *r  ■  0  and  h  S  0,  since 

We  may  conclude  that  unless  0  -   0,  <f>  ■  0  and  h  S.  0,  v  (x)  Sr  0  on 
(a,b);  i.e., 

By  an  analogous  argument,  the  function  v_  S  z.  -  u  satisfies 

and 

If  either  h  5p  0  or  not  both  &  and  <P>   are  zero,  then  we  may  conclude 

-?2(K)  ±  Hit)  £*   *t(i) 

on  [a,b]. 

Remark:   If  the  special  case  h  *  0,  &  -    <p  =0  holds,  then  from 
the  proof  of  Theorem  10  we  have  that  the  functions  v.  and  v  are  both 
identically  positive  constants,  say 

tfft)  s  M  f>  o     and  v±  Ms  J£z  >  o. 

Then 

i.e.,  we  have 

and  z.  and  z_  differ  by  a  positive  constant. 
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We  drop  the  restrictions  that  h  ^  0  and  that  6»  ,  <f>   are  in  [0,   /2], 
From  the  discussion  on  the  background  of  the  generalized  Sturm-Liouville 
boundary  conditions  we  may  have  just  as  easily  taken  €>  ,  ^»  in  (-fi"  fl/^/l] 
with  no  loss  of  generality. 

Theorem  11:   Suppose  that  u  is  a  solution  of  (7)  satisfying  boundary 
condition  (8)  with  -  1K/2  <     &     ^^/2,   -^/2  <  <j£  &  1^/2.   Let 
the  functions  z   and  z   satisfy  the  inequalities  (9),  (10)  and  (11),  (12) 
respectively  as  in  Theorem  10.   Then  the  bounds 

hold  in  (a,b)  if  and  only  if  there  exists  a  positive  function  w  on  [a,b] 
satisfying 

d3)      (lt+4t)lntj  £  0     0f     ^>V 
(14)      i4((h)aa><fi    +  cJtb)  Bfk<f   ±0 

in  such  a  way  that  not  all  the  inequalities  in  (13)  and  (14)  are  equalities, 

Proof:   Suppose  there  exists  a  positive  function  w  on  [a,b]  which 
satisfies  the  inequalities  (13)  and  (14).   Then  define  the  function  v  = 
u/w.   As  in  Theorem  4 

(W)LiQ=  ( 'L-H,)  IwJ 

which  we  may  rewrite  as 

(is)  (Z+#)[r]*  r"+  Qr1  +  Mr  »  & 

where  G  =  (2w"/w)  +  g  and  H  =  l/w(L+h)  w.   Note  that  H,  the  coefficient 
of  v,  is  non-positive.   The  functions  v  satisfies  the  boundary  condi- 
tions : 
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-jy  fy)  m )  +  Y(*J  toh)\  Cos  &    +  irti)  W(K )  fm  b  -  ^ 
\y*(h)  Uj/b)  t  v(b)  tv/^cJ^    +-  ?/«?)  ^W  j£i  ^  *  #, 


which  we  may  rearrange  to  form 


(16) 


Dividing  through  by  w(a)cos  (9-  and  w(b)cos<*»  ,  then  multiplying  by  the 
non-negative  quantities  cos  ii~r  and  cos  <p  in  the  first  and  second  condi- 
tions of  (16)  respectively,  we  have 

,                       [code  (-Lo'ra)  <xxe  +  t*/*J  &»&  )1       %  <3w & 
-vwcej-g  +  M)\ /=    — 

(17) 

r"*i  cts4  +  vt*)  \  — - — /  -   — — 

where  the  angles  &    and  ^  are  defined  by 


$m  a  •=. 


S^K)^" 


Mb)    O0Sf 
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or   equivalently 

to  (a)  Cose 

/1Q\ 

■f<Wd>    s - 

^  -TtfCh)  COS  4 

By   the   hypotheses   on     &  ,  *jP    and  w,    both   tan  d£    and   tan  <^>    are   non- 
negative.      By    (14)   we  may  always   take  3    an^  V      such   that 

Note:      If  9-    =  V/2,    take   jp   =  IV  /2 ;    if  <^  =  ft*/2   take   £    =  "ft* /2. 
Hence  we  have 

(19)  __ 

with    e"  »  ?    in   [0,    'ft'  /2]    and  ^j    «    ^  ^     cos  5)  /w(a)cos©    and      ^    = 
^r  cos   $2>  ^    /w(b)cos  <p 

If  the  functions  z  and  z   satisfy  the  conditions  (9),  (10)  and 
(11),  (12)  respectively;  then  the  functions  z  /w  and  z  /w  satisfy  the 
analogous  conditions  with  respect  to  equation  (15).   Hence  by  Theorem  10 
the  following  inequality  holds 

unless  £)  =  <f>    -   0  and  H  £  0.   From  the  proof  of  Theorem  10  this  situa- 
tion holds  when  the  inequalities  (9),  (10)  and  (11),  (12)  are  equations. 
Hence,  if  there  is  a  positive  function  w(x)  satisfying  (13)  and  (14)  but 
such  that  not  all  are  equations,  we  may  conclude 

(20)  2jt)  *  U(Z)  ^(X) 
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Note:   If  w  satisfies  the  boundary  value  problem 

(l+h)M  ^0      0*   fab) 

then  let  &  be  any  real  number  and  u  a  solution  of  (7)  which  satisfies 
(8).   Then 

fWJfti+Sv]*  fat)!*]  *  $M)M  =o 


and 


Thus,  if  w  satisfies  (13)  and  (14)  as  equations,  then  any  multiple  of 
w,  when  added  to  a  solution  of  (7)  and  (8),  will  yield  another  distinct 
solution.   Hence,  if  at  least  one  solution  u  exists,  then  there  are 
many. 

Now  consider  the  boundary  value  problem 

(2D      (l^h)[w]  «o     fa    (a,h) 


with 
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By  applying  the  classical  existence  theorems,  let  w  be  a  solution  of  (21) 
satisfying  (22).   Noting  that  the  function  w  satisfies  analogous  condi- 
tions to  the  function  u  in  Theorem  10,  we  have 

^  (%)  ±    IA//X)  ±   -2,  (KJ 

for  appropriately  chosen  bounding  functions  z  (x)  and  z  (x)  .   Also 
z-(x)  5£  0   satisfies  condition  (11),  (12);  hence  the  function  w(x)  is 
non-negative  on  [a,b].   We  now  show  that  w(x)  is  actually  positive  on 
[a,b]. 

If  w  is  zero  at  an  interior  point,  then  w'  is  also  zero  there, 
otherwise  w  would  have  some  negative  values.   But  the  identically  zero 
function  satisfies  the  same  initial  conditions  as  w(x)  and  by  Theorem 
5  for  uniqueness  of  initial  value  problems,   w  SE  0  on  [a,b]  which  con- 
tradicts (22).  Hence  w  cannot  vanish  at  an  interior  point. 

If  w  vanishes  at  an  endpoint,  say  at  a,  then  from  the  first  condi- 
tion of  (22)  we  have 

a/  VSrj  cos  &   a  -  i  . 

Now  w  non-negative  implies  that  w'(a)  £  0,  a  contradiction  since  cos^ 
^    0.   Thus  w(a)  >     0.   Similarly  w(b)  >  0,  and  we  have  w  ">     0 
on  [ a , b ] . 

As  with  Theorem  1  and  Corollary  1,  it  is  often  possible  to  prove  a 
more  general  theorem  and  by  a  judicious  choice  of  parameters  obtain  a  re- 
lated, more  specific  result.   For  example,  if  h  «  0,  0  ^  &  —  «*/2, 
0  £  <f  £   "IK/2;  then  the  function  w  <£  1  satisfies  (13)  and  (14)  and 
Theorem  10  follows  immediately  from  Theorem  11. 

The  fact  that  the  function  w  does  not  appear  in  the  inequality 


*i(t)  £  -v/k)  £  a,  ft) 
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prompts  us  to  try  to  obtain  a  necessary  and  sufficient  condition  that 
guarantees  that  the  functions  7.     and  z   form  upper  and  lower  bounds  for 
u  independent  of  a  function  w. 

Theorem  12:   Suppose  that  the  functions  z  (x)  and  z  (x)  satisfy  in- 
equalities (9),  (10)  and  (11),  (12)  respectively,  in  such  a  way  that 
equality  does  not  hold  in  all  the  conditions.   Let  g(x)  be  bounded  be- 
low on  each  interval  [a,c]  and  bounded  above  on  each  interval  [c,b] 
with  a  ^  c  S    b.   Let  u(x)  be  a  solution  of  (7),  (8).   Then  the  bounds 

(23)         $t\t)    ±   H(l)    £  TL,tl) 

hold  if  any  only  if  z_(x)  ^  z  (x)  for  a  <*  x  /  b. 

Proof:   Trivially,  if  (23)  holds,  then  z  (x)  £  z  (x) .  Hence  assume 

z      -  z   is  non-negative  and  we  must  show  that  (23)  holds.   Put 
2  (%)    =  ■£,  (jL)   -  ^  (I) 

If  q  ^  0  on  [a,b],  then  q  satisfies  the  conditions  for  the  positive 
function  w  in  Theorem  11;  i.e., 

(24)  (Lih)l^  =  fahjfhj-  &+ti)J*il  ^o  Oh  fah) 

since  z   and  z   satisfy  (9)  and  (11)  respectively. 

(25)  .  -, 

since  z   and  z   satisfy  (10)  and  (12)  respectively.   Hence  we  may  con- 
clude by  Theorem  11  that 

Thus,  we  need  only  consider  the  case  where  q  has  a  zero  in  [a,b].   We 
will  show  that  q  cannot  vanish  at  an  interior  point  and  that  if  q 
vanishes  at  one  or  both  endpoints  then  we  may  find  suitable  positive 
functions  w(x)  such  that  Theorem  11  will  hold. 
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Suppose  that  q(c)  =  0  for  c  in  (a,b).   Then  q'(c)  =  0  since  q  ^.  0 
on  [a,b],  i.e.,  q  has  a  local  minimum  at  c.   We  may  conclude  that  q  £j  0 
on  [a,b]  by  the  uniqueness  of  solutions  to  the  initial  value  problem. 
But  then  equality  holds  in  (24)  and  (25),  a  contradiction. 

We  are  left  with  q  ^  0  interior  to  [a,b]  and  q  vanishing  at  an  end- 
point,  say  at  x  =  a.   Then  q'(a)  ^>     0  by  above.   Thus  &    ='0^/2,  other- 
wise the  first  condition  of  (25)  would  be  violated.   Similarly  if  q(b)  =  0, 
then  q'(b)<  0  and  <£>    »  IK  /2. 

We  now  consider  the  case  where  q(b)  =  0  and  q(a)  ^  0,  and  proceed 
to  find  a  positive  function  w  such  that  Theorem  11  will  hold.  From  the 
preceeding  paragraph  <f>    -    Ir   /2.   Consider  the  initial  value  problem 

C  (L+h)[Y3   =  0 

I  ?Ya)  »  Ccs &      ansl  Y ?*)  -  Si'k  & 

which  has  a  solution  by  classical  existence  Theorems  for  differential 
equations.   Let  r  be  the  desired  (unique)  solution.   We  note  that  for 
6  <    "iP/2,  r  is  (strictly)  positive  at  a,  and  if  Q    =  *V  /2,  r(a) 
■  0  and  r  is  positive  near  a. 

Now  define  the  function  v  =  r/q  for  a  <  x  <  b.   Also  note  that 

rfa)    v 
j to)  3w  &  -  j'faJCQjQ 

by  the  first  condition  of  (25).   Using  the  facts  that  (1+h)  [r]  =  0 
and  (L+h)  [q]    0,  we  may  conclude  that  v  satisfies 
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with  G  =  (2q*/q)  +  g  and  H  =  l/q(L+h)  [q]  ^  0,  as  in  the  proof  of 
Theorem  11.   The  coefficients  G  and  H  satisfy  the  restrictions  in  Theorems 
2  and  3  on  any  subinterval  of  the  form  [a,c].  Hence  either:  (i)  v(c) 

y    v(a)  and  v1  (c)  ^>  0,  or  (ii)  v(x)  S  v(a)  for  x  in  [a,c].   In  either 
case  r(x)  >  0  for  x  in  [a,b)  since  both  v  and  q  are  positive. 

If  v(x)  is  identically  the  positive  constant  r(a)  on  [a,b],  then  q 
is  proportional  to  r  and  hence  must  satisfy  the  same  equations.  Thus 

and 

implies  that 

-  <g'fa)  Cos  &  -h   ffa)  $iM  &  =r  #. 

By  hypothesis  on  z  and  z   (hence  q)  equality  does  not  hold  in  all 
conditions  (24)  and  (25).   We  must  then  have 

from  the  second  condition  of  (25)  and  that  <fr  ■  1h  /2.  Hence  q(b)>  0, 
a  contradiction.  Thus  we  are  left  with  case  (i),  i.e.,  for  some  number 
c  in  (a,b), 

irtc)  >  ytej  ±o    and   ir'fjc)  >o     jcr  t  -c  . 


Now  *'*  -  W 

Put  f    (x)  s  r'q  -  rq1.   Then  from  (24)  and  (26) 

(27)  -v(^h)Lfij  =  -?*f"-  T'S?''-  Yh%  -o 


and 


us)    sMJM  'jrUp*  +  jhr   =£>. 
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Solving  for  -qhr  in  (28)  and  substituting  into  (27)  we  see  that 
satisfies  the  differential  inequality  p  -L^P  -  ~^  ( L'+hjtfl  -0. 
On  the  interval  [csb]   g  £  M  by  hypothesis,  thus 


Solving  for   Jr   we  may  write 


m=  fie 


and    for  x  ■  c 
thus 

fit)  k.  frc)e 


In  particular 

_  M  (t-a) 

since  v(c)  )>   0.   i.e.,  n 

But  q(b)  =  0  and  q'(b)  £r  0,  therefore  r(b)  >  0.   Hence  we  have  r 
positive  on  (a,b]  and  satisfying 

with 

-yfa;  cos s>  ■/■  r/aj  sw  o   ~£> 

and  since  <f>   =  ^  /2 
Hence  the  function 

tort)  s  3  w  -^  >vfcj 
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satisfies  the  requirements  oi  Theorem  11,  and  we  have 

On  the  other  hand,  if  q(a)  ■  0  and  q(b)  ^>  0,  analogously  we 
may  show  that  the  solution  of 

5(b)  =  cos4    <wd  sYi/-  ~sm<4 

is  positive  on  [a,b),  so  that  the  function 

satisfies  Theorem  11. 

Finally  if  q(a)  =  q(b)  =  0,  we  have  r  )  0  on  (a,b]  and  s  ^  0 
on  [a,b)  with  r(a)  =  s(b)  =  0;  hence 

^Wb  7YM)  +S(X) 

satisfies  Theorem  11.   This  concludes  the  proof  of  Theorem  12. 
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VI.   APPLICATIONS  OF  APPROXIMATION  FUNCTIONS 

Through  the  preceeding  sections  we  have  seen  that  the  results  have 
been  essentially  theoretical  and  that  examples  have  not  always  been 
fruitful  or  plentiful.   Sections  1-4  present  one -dimensional  maximum 
principles  and  their  use  in  some  differential  equations  theory  while 
Section  5  actually  gives  some  practical  techniques  to  apply  in  boundary 
value  problems.   The  problems  we  are  considering  are  second-order  dif- 
ferential equations  with  variable  coefficients  which  in  general  must  be 
solved  by  series  or  approximation  methods. 

There  does  not  seem  to  be  a  "best"  form  for  an  approximation  func- 
tion since,  in  general,  each  boundary  value  problem  needs  individual  at- 
tention.  However,  there  is  a  general  procedure  to  follow  to  determine 
bounds  and  behavior  patterns  of  the  solution  (under  certain  conditions) 
and  then  the  approximate  value  of  the  solution  at  any  point  in  the  inter- 
val under  consideration: 

1.  Find  the  straight  line  z(x)  which  satisfies  the  boundary 
conditions  as  equalities  and  apply  the  operator  (L+h)  to  it.   If  (L+h) 
[z]  is  of  one  sign  on  the  interval,  then  z  will  serve  as  either  an  upper 
or  lower  bound  as  determined  by  the  sign  of  (L+h)  [z],   (i.e.,  (L+h)  [z] 

&    0  implies  that  u(x)  £  z(x)  and  (L+h)  [z]  £   0  implies  that  z(x) 
^L  u(x)) . 

2.  If  the  constants  0     »  Q0   of  the  boundary  conditions  are 
of  the  same  sign,  then  the  approximation  function  z(x)  S  0  will  show 
that  the  solution  u(x)  is  either  non-negative  or  non-positive.   (i.e., 

&1    >     0,   'Jf  ^    0  implies  u  ^     0  and  ^     £    0,  /fr"  ^   0 
implies  u  £•  0) . 
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3.  If  the  general  behavior  of  u(x)  is  desired,  it  is  only 
necessary  to  bound  u  between  the  straight  line  from  (1.)  (assuming  this 
is  possible)  and  some  constant  function  which  satisfies  the  approxima- 
tion boundary  conditions  (9)  and  (10),  or  (11)  and  (12)  in  Section  V. 

4.  If  an  approximate  value  of  the  solution  u(x)  is  desired 

at  some  point  in  the  interval  then  find  a  polynomial  of  exponential  func- 
tion that  satisfies  the  approximation  boundary  conditions  mentioned  in 
(3.)  above.   If  a  more  exact  approximation  is  desired,  then  the  form  of 
the  approximation  function  is  more  complicated  and  the  difficulty  of  evalu- 
ating and  checking  the  restrictions  increases.   In  many  instances  it 
would  be  as  easy  to  find  an  approximate  or  series  solution  (if  initial 
conditions  are  known)  and  evaluate  it  numerically  on  a  computer.   The 
value  of  the  above  technique  is  that  the  behavior  of  the  solution  of  a 
complicated  second  order  differential  equation  can  be  determined  easily 
and  in  less  time  than  other  methods. 

The  following  examples  will  help  point  out  the  applications  of 
Section  5  and  Theorem  12  specifically.   The  choice  of  the  approximation 
functions  may  seem  ambiguous,  but  they  readily  become  more  understandable 
in  the  context  of  each  problem. 

Example  1:   Consider  the  boundary  value  problem 

(i)        K"  t  3V '-  TClt  *  O  OH     fatl) 

with 

(2)      -  v'to)  +  w(t>)  -  o      and       ?/ 0)  »  / . 

As  noted  in  (2.)  of  the  preceeding  discussion  the  solution  of  (1)  satis- 
fying (2)  is  non-negative.   Consider  the  straight  line  £i'X)  ~     %  ^^' ) 
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which  satisfies  (2).   Then 


on  (0,1) .   Hence  let  z(x)  ■  z  (x)  and  we  have 

We  note  that  y  S  1  satisfies  the  requirements  for  an  upper  bound,  i.e., 

O 
* 
(L+h)  [y]  £  0  and  y(l)  =  1,  -y1  (0)  +  y(0)  =  1  >  0.   We  have  quite 

simply  determined  the  general  behavior  of  u  in  (0,1),  i.e., 
(1/2) (x+1)  ^  u(x)  ^  1. 


##*3j 


k 


i/t)*  itx+i) 


iM 


Fig.  2 
To  determine  more  accurately  the  bounds  for  u  at  any  point  x,  we  choose 
a  more  complicated  function  for  z  ,  the  upper  bound.   A  natural  choice 
to  consider  would  be 


zx(x)  -  (1/4) (x+3) 


which  satisfies 


a(l)  -  1  and  -z'(0)  +  z(0)  =  3/4  >  0. 


But  z(x)  does  not  satisfy 


o)    (l+h)  [?.,-]  -  f[s-  *a+i)]  ±o 
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on  the  entire  interval  (0,1),   ba„  *c   may  define 

*  w  -  y 

where  (1/2) (  V  21-3)  is  a  zero  of  3  -  x(x+3)  =  0.   Then  z  is  continu- 
ous and  satisfies  the  conditions  for  an  upper  bound  and  we  have 

{  (k+i)   ±  lid)  &  i[\T\  +3] 

on   (0,(l/2)(    VTl   -   3)) and 

±(KH)     £   Ufa)       ±     \lpijf) 
on   ((1/2)    (      {Tl   -   3),1). 
Example  2;      Consider 

(4)  uM   -    (l-hc2)u  =  0   for  0  <    x  <    1 
with 

(5)  u(0)   =   1  and  u(l)    =  0. 

The  solution  u  must  be  non-negative.   Consider  the  straight  line 

zx(x)  =  1-x 
which  satisfies  (5).   Now 

*  0  —  <r/-txz)  o-x)  &  o 

since  0  <^  x  <£  1.   We  now  have  an  approximate  idea  of  the  shape  of  u. 
We  pick  for  z_  (the  lower  bound)  a  parabolic  function.   We  choose  z_ 
such  that  z_(0)  ■  1  and  z  (1)  =  0,  and  so  that  z   is  concave  upward. 
In  Example  1  we  showed  steps  1-3  of  the  general  method  and  here  we  show 
how  a  better  bound  can  be  found  by  choosing  a  more  complicated  approxi- 
mation function.   Let 

2 
z2(x)  -  (x-1) 
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which  satisfies  (5),   Then 


Hence 


*  Z-  (/+.**)/*- *JX  hoc*  /fjt) 

(£-0  *  ^  *M  £  /-X 


In  particular  x  -  1/2 

i  6  Mi)  &  %. 

Choosing  for  T_  the  expotential  function 
which  satisfies 


Now 


(M)fsfl*  **-  #***;*€ 


-x  3    2 

since  e    *  0  and  2  ^  x   -  x  .   Hence 


cf  *#-*>  ^  #/*;  ^  /-*, 


and  in  particular  at  x  =  1/2 


303  £  ufx)  £  .500, 


a  considerable  improvement. 
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Since  z,  (0)  =  Z-(u)  =  u(0),  we  may  bound  u'(0).   Now 

and  from  our  previous  results 

*•      2 
Example  3;  Consider  the  problem 

(6)        V0  -K*H    wp       -Jt^    0  <K  <  I 


with 


(7) 


-  V'fo)  t  h/c)  ^o    **d   7t(0  -  / 


Noting  the  similarity  of  boundary  conditions  (7)  with  boundary  condi- 
tions (2)  in  Example  1,  we  choose 

&t)m      i(K~+t) 
and  note  that 

Hence  z  (x)  =  (1/2) (x+1)  is  an  upper  bound  for  the  solution  u(x).   We 
choose  for  z   the  exponential  function 

which  satisfies  (7).   Then 

(i*+h)L*D  m    (t-t)*      ±o    #i  fat)  . 

Thus 

Xj-I 


f'  &  Ufa)  ±    i  (K+0 
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and  in  particular  at  x  =  1/2S 

0.6065  ^u(l/2)  ^  Oo7500. 
An  interesting  application  of  the  theory  developed  allows  us  to 
bound  solutions  to  the  differential  equation 

on  (a,b)  even  when  the  coefficients  g(x)  and  h(x)  are  piece-wise  continu- 
ous. 

Example  4:   Consider 

(8)  ^"t^*'-  h(t)K=o      Jer    #±*^l 

with 

(9)  y/cj^c      and    -ptO)  «•  / 

where 

f  1    ;    0  <x    ^  1/2 

h(x)    =  V 

(   2    ;    1/2  <    x     6.    1 

From  the  boundary  conditions  consider  the  approximation  function 

z(x)  -   x 
which  satisfies  (9).   Then 

Hence  0  ^  u(x)  ^  x9  since  u  must  be  non-negative  by  (9).   A  better 

lower  bound  may  be  determined  by  considering  the  parabola 

2 


z2(x)  -  x 


which  satisfies  (9)  and 
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Thus 


x  £  u(x)  £.  x. 


When  working  through  the  theoretical  background  for  this  technique, 
there  is  a  question  of  generalizing  the  approximation  functions  so  that 
the  "best"  bound  for  the  solution  u  can  be  obtained.   In  considering  this, 
the  following  discussion  will  point  out  some  problems  which  arise. 

We  will  choose  general  expressions  containing  parameters  for  the 
approximating  functions  z  and  z  and  attempt  to  pick  values  of  the 
parameters  to  make  z  and  t.~   satisfy  the  appropriate  conditions.  For 
example  choose  .  (jo—a  ) 


and 


where  A,B  and  C>(  are  constants  to  be  determined.   From  Section  V  we 
know  that  z  (x)  must  satisfy 

(l+tl)fa]    ±   -ftP)        on  (a.b) 


and 


Then 
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Now  choose   °(      so  large  that 

where  f,g  and  h  are  assumed  to  be  bounded  on  [a,b].   To  make  z  (x)  an 
upper  bound ,  we  must  have 

or 

A  [2  h/x)  -  pftfc  4  %(&+  h  m)  e         J  ±  -fit)  i 

This  is  satisfied  if  we  let 

^*  torn    l-2hit)+  (sz-*t<}M  +  ttfr))e  I 

and  then  if  we  choose  A  such  that 


By  its  definition  k  ^  0.   We  have  now  satisfied  the  condition 

f/iMJ  fa 7  -  J**)       ^v    #  c  x  <  t>. 

We  must  also  satisfy 

zff*)=  A 12-  e  J  *  A  >T>  j 

^/^-    A/  2-  e  J  £-*i .  . 

Li  <x(b~A) 
Let  c  ■  2  -  e       ,  and  choose  A  such  that 

i.e.,  A  is  the  largest  of  the  four  numbers  in  the  brackets. 
Then 

u(x)  -£  z, (x). 
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Using  the  same  "*(   in  z  (x),  we  choose  B  as  the  smallest  of  the  four 
numbers 


B  =  min  \%   >  S    W^K  |— f*«7  )  £1     . 

on  [a,b]. 

Example  5;   We  solve  the  following  problem  in  two  ways.  Consider 

u"  -  xu  -  0  for  0  <  x  <  1 
with 

u(0)  -  0  and  u(l)  -  I. 
Let  z  (x)  ■  x  which  satisfies  the  boundary  conditions,  and 

(L+h)    [z  ]  =  -x2  ^  0  on  (0,1). 
Hence 

0  £?   u(x)  £   x. 
For  z_  we  choose 

which  satisfies  the  boundary  conditions,  and  B  ^  0.   Then 

UJhjUil-  ■*X+8f*+*b*-K>2k  0         on  (0,x) 

for  B  ^  1/2,  which  can  be  seen  by  considering  the  values  of  B  which 
make 


since 
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and  we  are  dealing  with  a  cubic  equation.   Hence  we  choose  B  =  1/2,  and 
we  have 

(1/2) (x+x2)  ^  u(x)  £.  X^ 
since  z  &   z   on  (0,1).   From  the  theory  developed  before  this  problem, 
choose  the  approximation  functions  as  follows: 


*,/*>■=  a[z-  £**] 


and 


Then  we  must  choose  ©^  so  large  that 


«*2-  #   >  O  on  (0,1) 


hence   choose     GK     =1.      Then 


k  =  min[2*  4-   0-t)t  '  J  =   1        (at  x  =  0)  , 


and 


-1 

c   -  2    -  e 


Hence 


A  =  Max 


l>ir~>  *^] 


and 


Hence 


B-min^— r,    ^0]     -        £>. 


o  £  tffr)  £  — iL~~  /*2-  e\ 
2e~t     L  J 


which  does  not  improve  on  our  previous  results. 

The  method  outlined  previously  and  illustrated  in  the  example  also 
leads  to  very  interesting  aspects  of  bounding  solutions  of  second-order 
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boundary  value  problems  with  general  boundary  conditions. 
We  may  find  the  straight  line 

do)     2(K)  *  <?.y  +  e^ 

which  satisfies  the  general  boundary  conditions: 

(id  < 

if  the  parameters  £)  and  <zL    satisfy  certain  restrictions. 
We  may  rewrite  (11)  as 

6,  (-  COS  9)    -a  (qa  +<$.)  si*  9  *  % 

or  equivalently 

which  has  a  non-trival  solution  for  c.  and  c.  if  the  determinant  of  co- 
efficients        f"  ^ 

z>  = 

is  non-zero  (assuming      *$*"  +  "Y^  =£      0) . 
Hence  if  £^  and   Cp    are  not  solutions  of 

(12)       5(n(&  +  4>)  ■=  (a-b)    Sihe  5^h  / 
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then  we  may  always  find  the  line  (10).   We  note  that  if   S     +  <f     ■  nfV 
and    &   =  nlh   ,     ^  =  nlV   for  (n  «  0,  1,...),  then  (12)  is 
satisfied. 

In  Example  1  we  used  a  polygonal  line  as  an  upper  bound.   This 
procedure  may  not  be  fruitful  in  the  non-homogeneous  case.   The  reason 
for  this  will  be  made  clear  later. 

Example  6;   Consider  the  non-homogeneous  boundary  value  problem 

%"  -  It-    Zt% - K^ 

with 

1(fC)  =  1t(j)   =  O    . 

To  determine  possible  concavity  of  the  solution,  consider 

&*J   -  1 
then 

hence  z  ="  1  is  an  upper  bound  and  the  solution  u(x)  is  concave  upward 

since 

2-V=  -  )  4  2:  <0  OV\         (Oji). 

To  determine  a  lower  bound,  consider 

2/X)   *  3 
then 

hence  z  S  -3  is  a  lower  bound  since  the  appropriate  endpoint  in- 
equalities hold. 

By  substitution  we  may  verify  that 

is  the  solution.   We  might  expect  then  that  the  polygonal  line 
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22 1%)      m 


-  a  ;  o  4  it  £  ty 


would  satisfy  the  theoretical  conditions  for  a  lower  bound  since 


Fig,  4 
Note  that  this  is  the  opposite  of  the  procedure  in  our  theory  where  we 
determine  bounding  functions  by  applying  the  restrictions  to  them 
rather  than  determining  whether  actual  bounds  of  the  solution  satisfy 
the  restrictions.   Recall  that  for  z   to  be  a  lower  bound,  then 

(i)  (L-Hi)    [z  J  ^  f(x)  -  2  +  x  -  x  on  (0,1) 
(ii)  z  (0)  !==   0  and  z  (1)  ^  0„ 

Applying  the  operator  (L+h)  to  z,_9 


hi'  Mj*-K±  £ 


2-//-X2"  <M     f0J/)) 
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From  our  theory  this  would  seem  to  imply  that  z      is  an  upper  bound  for 
u;  but  by  an  oversight  we  have  not  considered  the  fact  that  each  segment 
of  the  polygonal  line  must  satisfy  the  restrictions  individually. 
Hence  for  z  (x)  -   -x  [09l/4]  we  see  that 

with 

Thus  our  theory  does  not  apply.   We  find  similar  results  for  the  re- 
maining segments  of  z_. 

This  points  out  an  important  result  that  for  the  non-homogeneous 
case  the  theory  developed  may  not  apply  to  the  polygonal  line  approxi- 
mation functions. 

The  reason  for  this  is  clear  upon  consideration  of  Example  1,  where 
we  choose  the  straight  line 

&*)*  iu+3) 

as  an  upper  bound 9  and  find  that 

changes  sign  on  (0,1),  i.e.,  (L+h)  [z]  has  a  zero  in  (091)  and  by 
suitably  changing  z(x)  we  can  satisfy  the  restrictions 

(MIL?!  -O     6M    fa/). 

But  for  the  non-homogeneous  case,  we  must  satisfy 

and  we  may  not  be  able  to  suitably  change  our  approximation  function 
to  satisfy  the  conditions., 

This  points  out  the  importance  of  parametric  approximation  func- 
tions and  being  able  to  choose  the  parameters  in  such  a  manner  as  to 
satisfy  the  restrictions. 
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Consider  then 


*tt)=  a[2-  e      J 
^,W  *  6  jfe  -  ^J 

and  choose  A,  B  and  o(  according  to  the  theory  developed  preceeding 

Example  5  so  that  z..  and  z?  are  upper  and  lower  bounds,  respectively. 

2 

Then  <o(     must  satisfy  ^  '     -    1    ^     0;  hence9  choose  ^  ■  2.   Then 

choose 

£  -  Af//?  /  2  +-  f  £    J  =  2  *  Se       *i '   K  «  ± . 

We  must  now  determine 


on 

and 


-  1 


Since  <J   ■  0,  we  do  not  need  c  ~  2  -  e   since   /q  /c  =  0.   Thus 


choose  A  and  B  so  that 


and 


L  z  +  3*  J 


-g-    nun    I  Oj  0 ,   _™-~~     0 


Thus 
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We  may  graph  the  above  by  approximating  the  function  z„(x),  and  we 
have 


-*/*) 


AM 


Fig.  5 

We  now  have  the  option  of  attempting  to  change  the  parameters  A, 
B  and  ©I  to  obtain  better  bounding  functions  z   and  z   or  keeping  the 
existing  functions  and  tolerating  the  possible  large  error.    In  either 
case  we  are  sacrificing  time  or  exactness  and  are  defeating  our  pur- 
pose. 

The  technique  used  in  the  preceeding  examples  has  indicated  the 
possibility  of  approximating  solutions  u(x)  by  multiplying  the  straight 
line  s(x)  found  above  by  "warping""  functions  £  (x)  and   ^  (x)  so  that 

This  technique  should  not  only  help  determine  the  general  behavior  of 
the  solution  but  yield  usable  numerical  approximations. 
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